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Non-specular reflections in a macroscopic system with wave-particle 
duality: spiral waves in bounded media. 
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(Dated: 3 April 2013) 

Spiral waves in excitable media possess both wave-like and particle-like properties. When resonantly forced (forced at 
the spiral rotation frequency) spiral cores travel along straight trajectories, but may reflect from medium boundaries. 
Here, numerical simulations are used to study reflections from two types of boundaries. The first is a no-flux boundary 
which waves cannot cross, while the second is a step change in the medium excitability which waves do cross. Both 
small-core and large-core spirals are investigated. The predominant feature in all cases is that the reflected angle varies 
very little with incident angle for large ranges of incident angles. Comparisons are made to the theory of Biktashev 
and Holden. Large-core spirals exhibit other phenomena such as binding to boundaries. The dynamics of multiple 
reflections is briefly considered. 



Wave-particle duality is typically associated with quantum 
mechanical systems. However, in recent years it has been 
observed that some macroscopic systems commonly stud- 
ied in the context of pattern formation also exhibit wave- 
particle duality. Two systems in particular have attracted 
considerable attention in this regard: drops bouncing on 
the surface of a vibrated liquid layei^ ^ and waves in 
chemical media^ ^^. The second case is the focus of this pa- 
per. We explore non-specular reflections associated with 
spiral waves in excitable media - reflections not of the 
waves themselves, but of the particle-like trajectories tied 
to these waves. 



I. INTRODUCTION 

Rotating spirals are a pervasive feature of two-dimensional 
excitable media, such as the Belousov-Zhabotinsky reac- 
^JQj^i3 16 pigure[l|a) illustrates a spiral wave from a standard 
model of excitable media discussed below. The wave char- 
acter of the system is evident. As the spiral rotates, a peri- 
odic train of excitation is generated which propagates outward 
from the center, or core, of the spiral. Much of the historical 
study of excitable media has focused on the wave character of 
the problem, as illustrated by efforts to determine the selection 
of the spiral shape and rotation frequency ^^tH] 

However, it is now understood that these spiral waves also 
have particle-like properties. This was first brought to the 
forefront by Biktasheva and Biktashe\^ and has been de- 
veloped in more recent years^ ^^ ^^^^^^. One of the more 
stri king il lustrations of a particle-like property is resonant 
JnypEElEi] shown in Fig. [T[b)-(d). Resonant drift can oc- 
cur spontaneously through instability, or due to spatial inho- 
mogeneity, or as here, by means of resonant parametric forc- 
ing (periodically varying the medium parameters in resonance 
with the spiral rotation frequency). As is seen, the core of the 



spiral drifts along a straight line. The speed is dictated by the 
forcing amplitude while the direction is set by the phase of the 
forcing, or equivalently the initial spiral orientation. 
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FIG. 1. Illustration of resonant drift and reflection for spiral waves in 
excitable media, (a) Periodically rotating spiral wave in the unforced 
regime. The wave rotates around a fixed core and the path of the 
spiral tip (white) is a circle, (b) Resonant drift. The medium is para- 
metrically forced at the spiral rotation frequency. The core moves 
along a straight path and the spiral tip traces out a cycloid (white), 
(c) Reflection of drifting spiral from a no-flux boundary. The incident 
and reflected angles, Oi and Or, are indicated, (d) Path of a drifting 
spiral in a square box. The underlying spiral wave at one instant in 
time is shown faintly. In all cases the plotted fields are the excitation 
variable, u, of the reaction diffusion model. Details are given later in 
the text. 

The trajectories of drifting spirals are unaffected by the do- 
main boundaries (or other spirals should they be present) ex- 
cept on close approach, where often the result is a reflection 
of the drifting core^^ ^^ ^^, as illustrated in Fig. [l[c). Reflec- 
tions are not specular - the reflected angle Or is not in general 
equal to the incident angle Oi. When placed in a square box. 



the drift trajectory typically will ricochet off each boundary in 
such a way to eventually be attracted to a unique square path 
where Oi -\- Or = 90°, as shown in figure [TJd). (This is the 
more common case, but others are considered herein.) 

The primary goal of this paper is firstly to determine ac- 
curately, through numerical simulations, the relationship be- 
tween reflected and incident angles for some representative 
cases of spiral waves in excitable media, and secondly to ex- 
plore the qualitative features of reflections in excitable me- 
dia, particularly multiple reflections in square domains. While 
the numerical and theoretical study of reflecting trajectories 
was undertaken by Biktashev and Holden many years agcP^^^Zl^ 
much more extensive results are now possible and desirable, 
especially since phenomena strikingly similar to that seen 
in Figs. [TJc) and (d) have been observed in other macro- 
sc opic system s with both wave-like and particle-like proper- 
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II. MODEL AND METHODS 



natural, unforced, spiral frequency. However, due to nonlin- 
earity, there is a change in the spiral rotation frequency under 
forcing and so ujf must be adjusted with A to produce reso- 
nant drift along a straight line. 

We have studied reflections in two situations. The first is 
reflection from a no-flux boundary. This type of boundary 
condition corresponds to the wall of a container containing the 
medium. We set the reflection boundary to be at x = and 
impose a homogeneous Neumann boundary condition there: 



du 
dx 



(0,y)=0. 



(4) 



Since there is no diffusion of the slow variable, no boundary 
condition is required on v. The medium does not exist for 
X <0. 

The second situation we have studied is reflections from a 
step change in excitability across a line within the medium. 
We locate step change on the line x = 0. We vary the thresh- 
old for excitation across this line by having the parameter b 
vary according to 



Our study is based on the standard Barkley model describ- 
ing a generic excitable medium^^. In the simplest form the 
model is given by the reaction-diffusion equations 



b{x,y) 



ifx> 
Ah ifx<0 



(5) 



—- = V^u + -u{l -u)[u 

at e \ a 

dv 



(1) 

(2) 



where u{x^y) is the excitation field (plotted in Fig. [T]) and 
v{x,y) is the recovery field; a, 6, and e are parameters. The 
parameters a and b collectively control the threshold for and 
duration of excitation while the parameter e controls the ex- 
citability of the medium by setting the fast timescale of exci- 
tation relative to the timescale of recovery. 

We consider two parameter regimes - known commonly as 
the small-core and large-core regimes. The small-core case is 
shown in Fig. [T] As the name implies, the core region of the 
spiral, where the medium remains unexcited over one rotation 
period, is small. This is the more generic case for the Barkley 
model and similar models and occupies a relatively large re- 
gion of parameter space in which waves rotate periodically. 
Small-core spirals are found in the lower right part of the stan- 
dard two-parameter phase diagram for the Barkley model (see 
figure 4 of Ref. 36 ). Large-core spirals rotate around relatively 
large regions [see Fig.[8ja) discussed below]. Such spirals oc- 
cur in a narrow region of parameter space^^ ^^ near the bound- 
ary for propagation failure. The core size diverges to infinity 
near propagation failure. 

Parametric forcing is introduced through periodic variation 
in the excitability. Specifically we vary e according to 



e{t) = eo (1 + A sm{ujft + (/>)), 



(3) 



where A and ujf are the forcing amplitude and frequency. The 
phase (j) is used to control the direction of resonant drift. The 
forcing frequency producing resonant drift will be close to the 



Unlike for the no-flux boundary, in this case waves may cross 
the line x = and so there is no boundary to wave propaga- 
tion. Nevertheless, drifting spirals may reflect from this step 
change in the medium and we refer to this a step boundary. 

The numerical methods for solving the reaction-diffusion 
equations are standard and are covered elsewhere^^ ^^. Some 
relevant computational details particular to this study of spiral 
reflections are as follows. A converged spiral for the unforced 
system is used as the initial condition. Simulations are started 
with parametric forcing and the spiral drifts in a particular di- 
rection dictated by the phase (j) in Eq. ([3j. The position of the 
spiral tip is sampled once per forcing period and from this the 
direction of drift, i.e. the incident angle Oi, is determined by 
a least- squares fit over an appropriate range of drift (after the 
initial spiral has equilibrated to a state of constant drift, both 
in speed and direction, but before the spiral core encounters a 
boundary). Likewise, from a fit to the sampled tip path after 
the interaction with the boundary, we determine the reflected 
angle Or. By varying (j) we are able to scan over incident an- 
gles. 

The simulations are carried out in a large rectangular do- 
main with no-flux boundary conditions on all sides. For re- 
flections from a Neumann boundary (|4]), we simply direct 
waves to the computational domain boundary corresponding 
to X = 0. We also study reflections more globally from all 
sides of a square domain with Neumann boundary conditions, 
such as in Fig. [ijd). In the study of reflections from the step 
boundary ([5]), the computational domain extends past the step 
change in parameter. We have run cases with the left com- 
putational boundary both at x = —7.5 and x = —lb and 
these are sufficiently far from x = that trajectory reflection 
is not affected by the computational domain boundary. The 
dimensions of the rectangular computational domain are var- 
ied depending on the angle of incidence. For Oi c^ ±90° we 



require a long domain in the ^/-direction, whereas for Oi ^ 0° 
a much smaller domain may be used. In all cases we use a 
grid spacing of h = 1/4. The time step is varied to evenly 
divide the forcing period, but At c^ 0.019 is typical. Except 
where stated otherwise, the model parameters for the small- 
core case are: a = 0.8, b = 0.05, and eo = 0.02. For the 
large-core case they are a = 0.6, b = 0.07, and eo = 0.02. 
For the step boundary bo = 0.05 and Ab = 0.025. Differ- 
ent values of the forcing amplitude and period, A and Uf, are 
considered. Given the desire to measure incident and reflected 
angles precisely, we have required drift be along straight lines 
to high precision and in turn this has required high accuracy 
in the imposed forcing amplitude and period. In the appendix 
we report the exact values for the forcing parameters used in 
the quantitative incidence-reflection studies. 



III. RESULTS 

Before presenting results from our study of reflections, it 
is important to be precise about the meaning of incident and 
reflected angles. As is standard, angles are measured with 
respect to the boundary normal. This is illustrated in Fig.[TJc). 
What needs to be stressed here is that spirals have a chirality 
- right or left handedness - and this implies that we need to 
work with angles potentially in the range [—90°, 90°], rather 
than simply [0°,90°]. 

Specifically, we consider clockwise rotating spirals and de- 
fine Oi to be positive in the clockwise direction from the nor- 
mal. We define Or to be positive in the counterclockwise direc- 
tion from the normal. Both Oi and Or are positive in Fig.[TJc) 
and for specular reflections Or = Oi. 



A. Small-core case 

We begin with the small-core case already shown in Fig.[T] 
Figures [2] and [3] illustrate the typical behavior we find in re- 
flections from both types of boundaries. In both figures the 
upper plot shows measured reflected angle Or as a function of 
incident angle Oi over the full range of incident angles. The 
lower plots show representative trajectories for specific inci- 
dent angles indicated. Here and throughout, the no-flux nature 
of the Neumann boundary is indicated with shading (x = is 
the at the rightmost edge of the shading), while the step in ex- 
citability at a step boundary is indicated with sharp lines. All 
parameters are the same for the two cases; they differ only in 
the type of boundary that trajectories reflect from. 

The reflections are far from specular. This is particularly 
striking for 6^^ < where the incoming and outgoing trajecto- 
ries lie on the same side of the normal. The reflected angle is 
nearly constant, independent of the incident angle, except for 
incident angles close to Oi = 90° . There is a slight variation 
in the reflected angle, seen as undulation in the upper plots, 
but the amplitude of the variation is small. 

One can also observe in the lower plots that the point of 
closest approach is also essentially independent of incident 
angle, except close to Oi = 90° where the distance grows. 



(b) 




FIG. 2. Illustrative results for reflection from a no-flux boundary, i.e. 
Neumann boundary condition, (a) Reflected angle Or versus incident 
angle Oi. (b) Representative tip trajectories showing reflections at 
the incident angles marked with crosses in (a). The reflected angle 
is nearly constant for the full range of incident angles. The forcing 
amplitude is A = 0.072. 



Spiral trajectories come much closer to the step boundary than 
to the Neumann boundary. 

It is worth emphasizing that there is no effect of forcing 
phase in the results presented in Figs. |2] and |3] As the inci- 
dent angle is scanned, the phase of the spiral as it comes into 
interaction with the boundaries will be different for different 
incident angles. While this could have an effect on the re- 
flected angle, we have verified that there is no such effect for 
the small-core cases we have studied, except at large forcing 
amplitudes near where spirals annihilate at the boundary (dis- 
cussed later). 

While we have not conducted detailed studies at other pa- 
rameter values, we have explored the small-core region of pa- 
rameter space. Figure |4] shows representative results at dis- 
tant points within the small-core region. The figure indicates 
not only a qualitative robustness, but also a quantitative insen- 
sitivity to model parameter values throughout the small-core 
region. In each case the upper plot shows 6^^ ^ 0° while the 
lower plot shows Oi ^ —70°. The reflected angle varies by 
only a few degrees throughout all cases shown in the figure. 
Case (a)-(b) is close to the meander boundary while (c)-(d) 
is far from the meander boundary and corresponds to a very 
small core. Cases (e)-(h) are relatively large values of param- 
eters a and b, both with Neumann and step boundary condi- 
tions. 

In the step boundary case, there is also the effect of A 6 to 
consider. Across a number of representative incident angles. 
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FIG. 3. Illustrative results for reflection from a step boundary, i.e. a 
step change in the excitability of the medium, (a) Reflected angle Or 
versus incident angle Oi. (b) Representative tip trajectories showing 
reflections at the incident angles marked with crosses in (a). The 
reflected angle is nearly constant for the full range of incident angles. 
The forcing amplitude is A = 0.072. 





(h) 



FIG. 4. Illustration of the insensitivity of reflections throughout 
the small-core region of parameter space. Upper plots show Oi ^ 
0°, including faint visualization of the i^-field at a particular time 
instance, while lower plots show Oi ^ —70°. The reflected angle 
is nearly constant independently of incident angle, parameter values, 
and boundary type. Model parameters span a substantial range of 
the non-meandering small-core region: in (a) and (b) a = 0.7, 6 = 
0.01; in (c) and (d) a = 0.95, h = 0.01; in (e)-(h) a = 0.95, b = 
0.08. Cases (g) and (h) are step boundary, the others are all Neumann 
boundaries. A = 0.072 throughout. 
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we observed that as Ah is incremented from 0.025 up to 0.05, 
the closest approaches of the spiral tips occur further from the 
boundary. We also find a slight reduction in the angle of re- 
flection. Decrementing A6 has the opposite effect. However, 
if A 6 is too small then the repulsive effect at the boundary will 
be too small and spiral cores will cross the boundary. 

We have examined the effect of forcing amplitude A. Fig- 
ures [5] and [6] show reflected angle as a function of incident 
angle for various values A as indicated. There is a decrease 
in the reflected angle with increasing forcing amplitude, or 
equivalently increasing drift speed. Generally there is also 
an increase in the oscillations seen in the dependence of re- 
flected angle on incident angle. The solid curves are from the 



Biktashev-Holden theory discussed in Sec.|IVA 



For sufficiently large forcing amplitudes small-core spirals 
may be annihilated as they drift into Neumann boundaries. In 
such cases no reflection occurs. We have not investigated this 
in detail as it is outside the main focus of our study on reflec- 
tions. Nevertheless, we have examined the effect of increasing 
the forcing amplitude through the point of annihilation for the 
case of a fixed incident angle ^^ = 0°. The results are summa- 
rized in figure [7] The reflected angle reaches a minimum for 
A c^ 0.11, and thereafter increases slightly, but does not vary 
by more than 4° up to the amplitude where annihilation oc- 
curs, A ^ 0.225, as indicated in figure |7ja). The forcing am- 
plitude at which annihilation first occurs is rather large in that 
it corresponds to displacing the spiral considerably more than 
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FIG. 5. Effect of forcing amplitude on reflection of small-core spi- 
rals for the case of a Neumann boundary. Points are measured re- 
flected angle as function of incident angle at forcing amplitudes A 
indicated. Solid curves are from Biktashev-Holden theory discussed 
inSec. HVAl 



one unforced core diameter per forcing period. Figure |7jb) 
shows tip trajectories on either side of the amplitude where 
annihilation occurs, while figure [tJc) shows annihilation at 
much larger forcing amplitude. We note that the exact ampli- 
tude at which annihilation first occurs depends slightly on the 
rotational phase of the spiral as it approaches the boundary. 
(Annihilation first occurs in the range 0.22 ^ A ^ 0.23 de- 
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FIG. 6. Effect of forcing amplitude on reflection of small-core spi- 
rals for the case of a step boundary. Points are measured reflected 
angle as function of incident angle at forcing amplitudes A indi- 
cated. Solid curves are from Biktashev-Holden theory discussed in 
Sec. HVAl 
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FIG. 7. (a) Reflected angle as a function of forcing amplitude A up 
to the point of annihilation at a Neumann boundary for small-core 
spirals. The incident angle is fixed at ^^ = 0°. (b) Tip trajectories a 
little below {A = 0.215) and a little above {A = 0.235) the forcing 
amplitude resulting in annihilation of the spiral at the boundary, (c) 
Tip trajectory at A = 0.05 showing annihilation at very large forcing 
amplitude. 



pending on phase.) Likewise, the spiral phase can affect the 
reflected angle by nearly 1° for A ^ 0.16. The influence of 
phase is nevertheless small for the small-core spirals. It is, 
however, more pronounced in the large-core case which we 
shall now discuss. 



B. Large-core case 

We now turn to the case where unforced spirals rotate 
around a relatively large core region of unexcited medium. 
This case is illustrated in Fig. [Sja) where a rotating spiral 
wave and corresponding tip trajectory are shown in a region 
of space the same size as in Fig. [T] The larger tip orbit and 
unexcited core, as well as the longer spiral wavelength, in 
comparison with those of Fig. [TJa) are clearly evident. While 
such spirals occupy a relatively narrow region of parameter 




FIG. 8. Illustration of a large-core spiral wave, (a) A portion of 
a rotating spiral and corresponding tip trajectory in a square region 
40 X 40 space units [same size as Fig.[T|a)]. (b) Resonant forcing 
and reflection for a large-core spiral shown in comparison to that of a 
small-core spiral. The forcing amplitude is A = 0.05 in both cases. 



space, they are nevertheless of some interest because asymp- 
totic treatments have some success in this region^^ ^^ and be- 
cause this is nearly the same region of parameter space where 
wave-segments studies are performed^ \^m ^ 

Figure [Sjb) shows a typical case of non- specular reflec- 
tion for a large-core spiral compared with a small-core spiral 
forced at the same amplitude. While many features are the 
same for the two cases, large-core spirals are found often to 
reflect at smaller dr and moreover, they can exhibit different 
qualitative phenomena. 

Figures [9] and [T0| summarize our findings for large-core 
spirals. Reflected angle as a function of incident angle for 
three forcing amplitudes is shown in Fig. |9] One sees the 
overall feature, as with the small-core case, that reflected an- 
gle is approximately constant over a large range of incident 
angles. This is particularly true of low-amplitude forcing, 
A = 0.022. However, there are also considerable differences 
with the small-core case. 

For large-core spirals the reflected angle increases with 
forcing amplitude. This is opposite to what is found for small- 
core spirals in Figs. |5]and|6] Moreover, the reflected angles 
are noticeably smaller than for the small-core case, as was al- 
ready observed in Fig.[8jb). 

We now focus in more detail on what happens in various 
circumstances. The left portion of Fig. [9] indicates the differ- 
ent dynamics we observe, depending on forcing amplitude, at 
large negative incident angles (Qi ^ —52°), and Figs. 10 a)- 
10 c) show representative trajectories with di ^ — 60 . At 
A — 0.05 trajectories glance off the boundary. That is, they 



remain close for short while before moving off with a well 
defined large negative reflected angle. The reflected angle is 
nearly constant at dr ^ —85° for incident angles di ^ —52°. 
At A = 0.036, di < -52°, trajectories become bound to the 



boundary and move parallel to it indefinitely. In Fig. [T0[c), 
with A = 0.022, one observes the trajectory moving along 
the boundary for a distance before abruptly leaving the bound- 
ary at a well-defined, relatively small positive reflected angle. 
This behavior is not restricted to ^^ ^ —52° and is observed 
until 6i ^ +20°. In fact, this type of reflection is also ob- 
served for the other two forcing amplitudes studied for Oi in a 
range above —52°. For A = 0.036 this occurs until Oi is ap- 
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FIG. 9. Summary of results for large-core spirals. Reflected angle is 
plotted versus incident angle for three forcing amplitudes as labeled. 
Neumann boundary conditions are used. For A — 0.05 spirals are 
frequently annihilated at the boundary, [Fig.[Tofd)], over the range 
of incident angles indicated. For A — 0.05 and Oi ^ —52° trajec- 
tories glance from the boundary [Fig. [Tofa)]. For A = 0.036 and 
Oi ^ —52° trajectories become bound to the boundary. [See text and 
Fig.fTofb).] Wiggles are the effect of incident phase. 




FIG. 10. Catalog of interesting trajectories for large-core spirals, (a)- 
(c) show impacts with Oi ^ —60° at different forcing amplitudes, (a) 
With A = 0.05 the trajectory glances from the boundary and moves 
off nearly parallel to it (^r ~ —85°). (b) With lower amplitude 
A = 0.036, the trajectory becomes bound to the boundary, (c) With 
yet lower amplitude A = 0.022, the trajectory hugs the boundary for 
a while then leaves abruptly at an oblique angle {Or ~ 23°). In (d) 
the incoming spiral with large forcing, A = 0.05, is annihilated at 
the boundary. In (e) the effect of phase is seen with two approaching 
trajectories shifted by half a core diameter. Otherwise the conditions 
are identical, A — 0.05. The resulting reflected angles differ slightly. 



proximately —15°, while for A = 0.05 this is seen only until 
i9i is about -45°. 

At the higher forcing amplitudes, as indicated for the case 
A = 0.05 in Fig. |9] large-core spirals are frequently annihi- 
lated when they come into contact with the boundary. Fig- 
ure [TOjd) shows a typical example. Whether or not a spiral 
is annihilated depends very much on the spiral phase on close 
approach to the boundary. The points shown in Fig. [9] with 
A = 0.05 are those where the trajectory reflected; the ab- 
sence of points indicates annihilation. However, these results 
are for spirals all initiated a certain distance from the bound- 
ary. Changing that distance would affect the spiral phase at 
close approach and hence a different set of points would be 
obtained. Nevertheless, the marked range of frequent annihi- 
lation is indicative of what occurs at this forcing amplitude. 

Finally, we address the wiggles in the reflected angle curves 
in Fig. [9j most evident at large forcing amplitudes. These 
wiggles are also due to the fact that the phase of spirals on 
close approach varies with incident angle. Figure [TOje) illus- 
trates how the reflected angle depends on phase by showing 
two trajectories shifted by half a core diameter. This shifts 
the spiral phase upon approach to the boundary and results 
in slightly different reflected angles. Rather than eliminat- 
ing these wiggles by averaging over various initial spiral dis- 
tances, we leave them in as an indication of the variability due 
to this effect. In general, reflections of large-core spirals are 
much more sensitive to phase than reflections of small-core 
spirals, and one should understand that the data in Fig. |9] will 
vary slightly if similar cases are run with spirals initiated at 
different distances from the boundary. 

While we have not studied the step boundary in detail for 
large-core spirals, we have carried out a cursory investigation 
for such a boundary with Ab = 0.035. With the exception that 
there is no annihilation at the step boundary, we observe qual- 
itatively similar behavior to that just presented for the Neu- 
mann case. Most notably we find both glancing and bound 
trajectories. 



IV. DISCUSSION 



A. Biktashev-Holden theory 



Many years ago Biktashev and Holderr^*^ carried out a 
study very similar in spirit to that presented here. More- 
over, they understood that a primary cause for the reflection 
from boundaries was the small changes in spiral rotation fre- 
quency occurring as spiral cores came into interaction with 
boundaries. Based on this they proposed an appealing simple 
model to describe spiral reflections. The model is based on 
the assumption that both the instantaneous drift speed normal 
to the boundary and spiral rotation frequency are affected by 
interactions with a boundary, with the interactions decreas- 
ing exponentially with distance from the boundary. While 
the actual interactions between spiral and boundaries are now 
known to be more complex (see below), it is worth investi- 
gating what these simple assumptions give. The beauty of the 
simple model is that it can be solved to obtain a relationship 



between reflected and incident angles, depending on only a 
single combination of phenomenological parameters. (They 
called this combination 9, but we shall call it p. They also 
used different definitions for incident and reflected angles.) 

The model naturally predicts large ranges of approximately 
constant reflected angle depending on the value of p. What 
is nice is that while fitting the individual phenomenological 
parameters in their model would be difficult, it is also unnec- 
essary. The value of p can be selected to match the plateau 
value of Or observed in numerical simulations. Then the en- 
tire relationship between Or and Oi from the theory is uniquely 
determined. 

Curves from the Biktashev-Holden theory are included in 
Figs. [5] and [6] While there are obvious limitations to the the- 
ory, it is nevertheless interesting to see that some of the fea- 
tures are reproduced just from simple considerations. The the- 
ory would be expected to work best where the drift speed is 
small: low amplitude forcing. For the large-core spirals the 
theory does not apply and so the corresponding curves are not 
shown in Fig. [9] 

We do not intend here to propose a more accurate theory 
for the reflection of drift trajectories. It is worth emphasiz- 
ing, however, that in recent years much has been understood 
about the interaction between spiral waves and various sym- 
metry breaking inhomogeneities^^^Hl jj^g j^^y ^q ^j^j^ under- 
standing has been response functions, which are adjoint fields 
corresponding to symmetry modes for spirals in homogeneous 
medi^^. In principle, one could now obtain rather accurate 
description of reflections using this approach, at least for cases 
in which the boundary could be treated as a weak perturbation 
of a homogeneous medium. This is beyond the scope of the 
present study and we leave this for future research. 



B. Multiple reflections 

As noted in the introduction, Fig.[TJd), when placed within 
a square domain the trajectory of a drifting spiral will typi- 
cally approach a square, reflecting from each domain bound- 
ary such that Oi ^ Or = 90° . The reasons for this are simple 
(see for example Prati, et al.), but a brief analysis is use- 
ful, particularly for understanding when square orbits become 
unstable. 
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FIG. 1 1 . Sketch showing the geometry of multiple reflections in 
a portion of a square box of normalised length, d^ is the location, 
relative to the length of a side, of the n*^ reflection. 



flections in the case where the reflected angle is larger than 
45°. In this case the path will necessarily strike consecu- 
tive sides of the domain. Consider first the path in terms of 
angles and let Of and 0'!^ denote, respectively, the n^^ inci- 
dent and reflected angles, starting from the initial reflection 
O^^O^. Trivially, the geometry of the square domain dictates 



that 0, 



n+l 



O'f = 90°. Then, if the trajectory approaches 



an attracting path with constant angles, lim^^oo ^^ = ^t^ 
lim^^oo ^? = 0*, it must be that this path satisfies 6>* + 6>* = 
90°. That is, it must be a square or a rectangle. Denot- 
ing the relationship between incident and reflected angle by 
Or = 0{Oi), then a necessary condition for the square path 
to be attracting is that |9'(6>*)| < 1. For the cases we have 
studied this is true since \Q\0*)\ c^ 0. 

Turning now to the points at which the path strikes the edge 
of the domain, we let d^ denote the position of the n^^ reflec- 
tion along a given side, relative to the length of a side. One 
can easily see from the geometry that (i^+^ = (1 — (i^) cot ^^. 
Now, since Or -^ 0*, the fixed point d* is given by d* = 
(1 - d*)cot(9;, or d* = 1/(1 + tan(9*). This corresponds 
to a square trajectory. For example, from Fig. [5] with a forc- 
ing amplitude A = 0.072 one can see that 0* will necessarily 
be about 74°, giving d* ^ 0.22. These are the values seen 
in the simulation in Fig.[TJd). A necessary condition for this 
fixed point to be stable is | cot 0^\ < 1. For small-core spirals 
6>* > 45°, so cotO* < 1, and hence their square paths are 
stable. 




FIG. 12. Examples of non-square paths for large-core spirals, (a) 
A = 0.022. The reflected angle is considerably smaller than 45° 
and the resulting trajectory bounces between opposite sides of the 
domain. The spiral is shown faintly at one time instance, (b) 
A = 0.034. The square trajectory is unstable. For the first circuit 
around the nearly square path the full tip trajectory is plotted. Subse- 
quently, for clarity only, the tip path sampled once per forcing period 
is shown. The trajectory collapses towards the diagonal. The final 
portion of the trajectory before the spiral approaches the corners is 
shown in bold. The spiral undergoes a complicated reflection from 
the corner (gray, dotted). 

While square trajectories occur for small-core spirals, for 
large-core spirals other trajectories are possible. Examples 



are shown in Fig. 12 These occur when the reflected angle 
is smaller than 45°. (It is possible that Or < 45° might occur 
for small-core spirals in some regimes, although we have not 
observed them.) When Or < 45° it is not necessarily the case 
that trajectories will strike consecutive sides of a square box. 



Figure 1 1 shows the geometry of a consecutive pair of re- 



This is seen in Fig. 12 a) where the spiral in reflects between 
opposite sides of the domain. The reflections satisfy Or = 



The more interesting case is when 0^ is only sHghtly less 



than 45° as is seen in Fig.[T2[b). The square trajectory is un- 
stable. While \Q'{0*)\ < 1 and the angles converge quickly to 
0"^ + 6>;? c^ 90°, the equation d^+i = (1 - d"") cot 0"^ exhibits 
growing period-two oscillations for cot ^^ slightly larger than 
1. Period-two oscillations in d'^ with 6^ ^ 6^ c^ 90° corre- 
spond to approximately rectangular trajectories that approach 
a diagonal. This ultimately leads the spiral into a corner of the 
domain where it may reflect in complicated manner. 

The analysis just presented should not be viewed as a model 
for spiral trajectories. Rather it just shows what global dy- 
namics can be deduced simply from a measured relationship 
between incident and reflected angles. Essentially this same 
analysis appears as part of a study of cavity solitions^^ which 
also undergo non-specular reflections from walls and hence 
exhibit square orbits similar to figure [TJd). Our simple anal- 
ysis should be contrasted with the situation for drops bounc- 
ing on the surface of an oscillating liquid, so called walkers. 
Here physical models of the liquid surface and drop bounces 
account for many varied features of the sy stenr^^^^ . The 
corresponding theory for spiral waves would be that based on 
response functions in which many details of drift trajectories 
can be predicted^ ^ ^^, although as of yet not reflections from 
boundaries. Memory effects are important for walkers be- 
cause bouncing drops interact with surface waves generated 
many oscillations in the past, and models necessarily take this 
into accounFE3H3 However, path memory is absent from spi- 
ral waves in excitable media and this constitutes a significant 
difference between the two systems. 



propagating wave segments studied by Showalter et al^^^^ 
or the drifting spirals in excitable media considered here. (We 
examined briefly small-core drift trajectories through a single 
slit, but did not observe diffraction-like behavior.) Neverthe- 
less, for the reflection problem, spiral trajectories, propagating 
wave segments, cavity solitons, and wal kers all share the fea- 
ture of non- specular reflectionP^^ESEH and as a result these 
systems can show similar dynamics wh en underg oing multi- 
ple reflections within a bounded regiorP^E^'^^'^. It will be 
of interest to make further quantitative comparisons between 
these different systems in the future and to explore theoretical 
basis of this behavior. 



APPENDIX: FORCING PARAMETER VALUES 

In the appendix we report the exact values for the forc- 
ing parameters used in the detailed quantitative incidence- 
reflection studies, since obtaining high-precision values for 
resonant drift can be time consuming. The values stated in the 
body of the paper are reported only to two significant figures. 



TABLE L Parameter values used to produce resonant (straight line) 
drift in the small-core case (a = 0.8, 6 = 0.05). 



A 



^f 



At 



0.044462 
0.071868 
0.102609 
0.196132 



1.82 
1.792 
1.75 
1.63 



0.0187625 
0.0188508 
0.0188968 
0.0188957 



C. Concluding remarks 

We have reported some quantitative and some qualitative 
features of resonant-drift trajectories in excitable media. The 
main message is that reflections are far from specular - the 
reflected angle generally depends only weakly on the inci- 
dent angle and typically is nearly constant over a substan- 
tial range of incident angles (pa rticul arly negative incident 
angles). Biktashev-Holden theory ^^^^ accounts for some of 
the observed feat ures, b ut a more detailed theory based on 
response functions^E^EH jg needed. We have seen that the be- 
havior of large-core spirals is more varied than that for small- 
core ones. Rather than simply reflecting from a boundary, 
large-core spirals may sometimes become bound to, or glance 
from, or be annihilated at a boundary, even at moderate forc- 
ing amplitudes. Finally we have considered what can occur as 
spirals undergo multiple reflections within a square domain, 
and in particular have shown that while small-core spirals are 
observed to meet the conditions of stable square trajectories, 
large-core spirals may fail to meet these conditions and exhibit 
more interesting dynamics. 

We motivated this study with a broader discussion of 
macroscopic systems with wave-particle duality. A large 
number of analogues to quantum mechanical systems have 
been reported for walkers on the surface of a vibrated liquid 
laye#^^. As far as we are aware, this is less the case for the 



TABLE XL Parameter values used to produce resonant (straight line) 
drift in the large-core case (a = 0.6, b = 0.07). 



A 



^f 



At 



0.022 
0.035863 
0.050144 



1.025 
1.003 
0.989 



0.0188035 
0.0187557 
0.0187961 
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